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ABSTRACT - 


An approximate boundary for a finite sample sequential decision 
process is presented (without derivation). By means of Computer Sim- 
lation, properties of the process are checked; the power of the test is 


thus determined for selected alternative hypotheses. 
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I. INTRODUCTION 


The main purpose of this paper is to study the properties of a 
mathematically derived sequential sampling procedure. The derivation 
involves the assumption that a certain parameter N, identified as a 
sample size, was "large." Our purpose was to make computer calculations 
that would reveal how the sampling procedure might behave for finite, 
or even rather small, N. Other calculations would also indicate the 
power characteristics of the test, that is, the probability of rejecting 
a null hypothesis. The null hypothesis, in this case, was that a partic- 
ular probability distribution describes the sample. By experimental 
sampling methods we attempted to estimate the power that the procedure 
possessed to detect departures of various sorts from the null hypothesis. 

The stochastic process used to obtain the sequential procedure might 
also be used for other purposes in operations research. For example, 
the sequential procedure might be used to describe arrivals in transitory 
Queueing situations. Some of these situations are briefly described at 


the end of this paper. 





IIT. A SEQUENTIAL SAMPLING PROBLEM 


Suppose, for example, that a supplier of transistors asserted that 
each of his transistors had an in-service lifetime T, where T was a 


random variable distributed exponentially; 


T-e M8” 420 
F(t) = 


O a 0a 

In order to test that assertion or hypothesis, a sample of N tran- 
Sistors was selected am placed in service simultaneously at time t = 0. 
On the average, the number of transistors that failed by time t >0 was 
NF(t), i.-e.- if A(t) denoted the number of transistors that failed by 
time t, then 

E fA, (t)} = NF(t) 
and 


Var fan (+ )} = NF(t) [2-F(+)} ’ 


Since the transistors were presumed to fail, or die, independently. 

If one waited until all N transistors had failed, the rate could be 
estimated and compared to 5° A hypothesis test can be conducted by 
waiting until r out of N transistors failed (r<N); this was shown by 
Epstein and Sobel [1]. 

The present approach to this problem was to establish a boundary, 
By(t), with the property that if A= A 52 i-e. if the supplier was en- 
tirely honest, then Ay(t) < By (t) for all t, with a known, specifiable, 
probability denoted by 1-6. For instance, suppose 6 = 0205. Thenneat 


A = Ao» the probability that the accumulated numbers of failures by 





time t would ever exceed the boundary is five per cent. Hence, if Ay 
oegeans below the boundary, A<A “9 with ninety-five per cent confidence. 

The detailed derivation of the boundary is given by Gaver, Lehoczky, 
and Perlas in an unpublished technical report that is in preparation. The 
result is quoted below. 

The Boundary 

If F(t) is sufficiently smooth (possesses a continuous density) and 
nate 

B(t).= NF(t) +V¥ {ar(t) +p Li-r(t}} , (2.1) 
then the probability that A(t) < BL (t) for all % is appromimavely 
‘eal 3; the approximation improving as N becomes large (N~@). Here, 
x and p are arbitrary positive Be ois. 

A graph of a possible failure arrival patter is a step function with 
unit increases at time points that are the order statistics of F(t): 
T1): T(2)s eee where T(1) is the smallest of the N failure times, T(2) 
is the next smallest, ... TK) is the kth smallest. Then, a particular 


sample may be represerted as follows: 


A(t) =O) for (0 = teeer 


(1) 
= | for aa) << T5) 
= 2 for T 5) = ts 13) (222) 
= k for ti) <t< Tet) 
= a] © tor TN) <b>, 


and graphically, we have the following picture. 








Ee 2) aS) ke 


Figure 1. 


iierrveure 1, a crossing of the boundary occurs at Toe We would conclude 
that the null hypothesis is not satisfied if such an event occurs. Notice 
that the parameters & and B may be selected at will. It is interesting 


a 0.05, 


to arrange that «Bp = k, a constant such that — 
for example, and then to attempt to select ( or B ) so that the 
boundary will be crossed soon if a non-null distribution of specific 
sort is to be detected yale This problem has not been extensively in- 
vestigated, but the expected time to hit the boundary has been tabulated 


for certain of the problems studied. 





TTI. USE OF THE COMPUTER SIMULATION 


Several facts about the boundary By (t) presented in equation (2.1) 
are worthy of notice. 

a. The prescribed probability of crossing the boundary, Ma ae 
is not exact for finite N. That is, use of this probability to assess 
the probability of crossing is an approximation which is, in fact, quite 
analogous to the use of the normal distribution as an approximation to 
the distribution of a sum of independent random variables. The mathema- 
tical derivation does not indicate whether the actual probability of 
crossing is, perhaps, consistently lower or higher than the actual true 
probability. More complex mathematical methods may provide useful in- 
formation of this type, but a straightforward approach is that of simu- 
lation or synthetic sampling. By vats means, it can be seen how increases 
in N improve the quality of our approximation, as is suggested by the 
mathematical derivation. 

b. The relationship of the true crossing probability to the 
approximate crossing probability, e~2X PB for given N, may possibly be 
influenced by the ee of « and Pp - Setting XB =k, where k is a 
constant delivering a desired probability of crossing (for 5%, k=1.h9), 
x ee hence PB) may be adjusted so as to obtain good agreement between 
the theory (oes =e ©*%P ) and the true probability for various values 
of N-- in particular for small N. In this paper, an attempt was made 
to study this problem by simulation. 

In order to study the properties of the boundary By(t), the following 


Simulation was carried out. 








A. SIMULATION LOGIC 
1. The null hypothesis distribution was chosen and called Fo(t)- 


The formla used to calculate the boundary values was as follows: 


By(t) = NEg(t) + VN Caro(t) +p [1-Fo(t]} - Baty 


An exponential distribution is the null hypothesis, as in the illustra- 


tive problem; hence, for the initial calculations, the equation 


By(t) = N(l-e) + VN {p+ (a -P) (2-e')} (3.2) 


é 


was used. 

2. A sampling distribution was selected. Initially, it was the 
null hypothesis distribution, so that the accuracy of the theoretical 
estimate of the probability of crossing the boundary could be checked. 
later, other distributions were tried in order to determine the power 
for departures from the exponential distribution. 


3. A sample of size N was taken and arranged in order: 


16 <t 


tis (2) oe vag 7) a ae (38 


This ordered sample was composed of the actual numbers obtained as 


realizations of the random variables. 


Ta, Teo) ST Seay: (rp) 


The details on how the random samples were conducted will be presented 
Shortly. 

kh. At-each failure time, say t(K)s starting with k = iy then 2, 3, 

ae An (tx) = k was compared with By (ty): ef kK< By(teyy); for k = 

1, 2, ..-, the process continued to k + 1 and was checked again. If at 


kx, the event kx 2 Bult cy) occurred for the first time, a boundary 


y 





crossing was recorded at the time t. & ty.3 the test was stopped and 
another simulation begun. 

5. Step four, above, was independently repeated NR times (NR means 
the number of replications). The probability of crossing ty was esti- 


mated by the fraction 


_ _number of crossings. 
ty = NR (3.5) 


Clearly, NRfy was a realization of a binomially distributed random 
variable. An approximate idea of the variability of the estimate of 


the probability of crossing is given by 


f_(1-f, ) 
6 : _ a i (3.6) 
NR NR 


An approximate confidence limit for the true probability Pe of crossing 


ne 
came from assuming that ty is approximately normal with variance Oa 


B. SIMULATION PROGRAM 

The simulation procedure just described was written in FORTRAN com- 
puter programming language and run on the IBM 360, Model 67 Computer at 
the Naval Postgraduate School. The program was written to facilitate 
changes of distributions, boundary parameters (« and Ble and numbers 
of replications (NR). It was also constructed so that the sampling 
process could be stopped after a fixed time period, i.e.- before all of 
the order statistics had been recorded. 

The pseudo random number generator employed to obtain realizations 
of uniform random numbers was a library subroutine called RANDOM [3 and lj; 
it is based on a Jehmer multiplicative scheme. The uniform realizations 
were then transformed into other random numbers by a probability integral 


transformation. 


10 





A program listing and sample output has been included at the end 
of this paper. In order to generate the sequence of order statistics 
(interpreted as failure or arrival times), the random sample of size N 
was first created. Then, a SORT routine was employed to put the times 
in order. A more efficient procedure would have been to generate the 


ordered increments sequentially; for example, the time between af ) and 


k 


ay in exponential sampling is also exponential with a mean propor- 


tional to One program has been written in this vein, but production 


pul 
N-k’ 
runs for the present investigation have largely been carried out by use 


of the SORT routine (5]. 


Vy 





IV. SIMULATION EXPERIMENTS 


Several simulation experiments have been carried out to develop an 


understanding of the boundary crossing process. 


A. EXPERIMENT 1 

This experiment explores the crossing probability under the null 
hypothesis ( A= ro) and under selected alternatives, for various sample 
sizes (N)s; and for different boundary parameters & and PB , but with 
A = Bp. The computed probability of a crossing fy was calculated by 
the fraction (3.5) mentioned previously. The average number of arrivals 
up to a crossing (A(t)) was calculated by summing the number of arrivals 
(Ayy(t )) that had occurred up to the time of crossing, over the two thou- 
sand replications and dividing by the number of angen (NC) recorded. 
The average crossing time (C(t)) was calculated by adding the arrival 
times, which generated crossings, then dividing that sum by the number 
of crossings. The results appear in Table I (theoretical crossing prob- 
ability of five per cent), and Table II (theoretical crossing probability 
of ten per cent). These Tables are based on two thousand replications 
(NR = 2000). 

When the boundary was not crossed for a sample size N, all values for 
that N appear as O. In such a situation, the zero (0) is intended as a 


meaningless place holder rather than the unique number. 
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TABLE I. 


5% Probability of a Crossing 
X= B= 1.22388 


EXPON ( A= .5) 


= 


NC A(t) C(t) fig 
3 Sree: . ee 0.1% 
2 5.00 300 0.10 
(i 6.2 28 0.35 
1 8.00 523 0.05 
0 0.00 000 0.00 
0 0.00 000 0.00 
0 0.00 .000 0.00 

A= 1) 

NC A(t) C(t) fy 

6 3.63 Zar 2.30 

69 5 Oy 391 3.45 

66 6.30 02 3.30 

70 OT 16 2.50 

s 15.19 aie BOS 

85 27.5h 50k W.25 
104 3 586 Cae 

-EXPON ( A= 1.25) 

NC A(t) C(t) a 
100 3.66 298 5.00 
11 oi 381 7.G5 
165 6.33 Hh 8.25 
197 7.8 28 9.8), 
302 15295 BS || 2 18.59 
6h) 28.64 52 32.19 
839 1.09 560 1.95 
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TABLE I. (continued) 
EXPON ( A= 1.5) 


NC A(t) c(t) f 


_ — = atl 
17h 3.68 . 303 8.70 
253 5.03 362 12.6), 
331 6.26 2395 16.5 
hoh 7.8 26 20.19 
896 11S ee 98 hh.79 

1461 A Ei 470 73.04 
1770 55 5S ~ 1g? 88.9 


NC 


262 
389 
525 
672 
1391 
1893 
1985 


NC 


360 


71 
Fed 
1736 
1985 
2000 


EXPON ( A= 1.75) 


EXPON ( A= 2) 


A(t) 


3.69 
1.98 
6.15 
fio 3 
13.95 
19.86 
23.65 


C(t) 


301 
349 
381 
18 
NST 
. 362 
Rein 


= 


oe) 
1s 
PAS PN 
33.59 
69.5) 
9) 6 
99.2h 


reens 
Ota 
37.04 
h6 uh 
86.79 
99 2h 


100.00 








= 


NC 


—> 


Se) On = NON 


NC 


127 
TS 


156 
160 
It 
190 


TASES Ti. 


10% Probability of a Crossing 
Pp = 1.07298 


ele 


EXPON 


> | 
oo" 
ct 
a 


OOO WVUIVIW 
@ Ga@ ©>—~) © 
® @oaonm ow 


EXPON ( 


EXPON ( A= 1.25) 


A) 


3.51 
vege 
5.98 
7.10 


15.1h 
26.86 
38.66 
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A= 


-5) 


C(t) 


- 300 
-467 
341 
. 4,00 
000 
000 
-000 


C(t) 


eo 
428 
nites 
les 
HOS 
522 
551 


c(t) 


-3L5 
381 


25 ; 


517 
ey 
513 


kt 


|2 


SSO Or. Oe 
S @@ iy = J 
OO OVILMVUIVNI 


a 


6.35 
75h 
(oe 
1-79 
8.00 
8.85 
ATIC 


i2ne5 
1.45 
Wes 76 
18.70 
Sash 
h7.2h 
58.39 





NC 


360 
160 
571 
67) 
1197 
1682 
1907 


NC 


498 
660 
823 
7 ite 
eels 
195K 
1998 


NC 


617 

868 
1083 
127 
1878 
1997 
2000 


TABLE II. 


EXPON ( A= 1.5) 


A(t) 


3 Sik 
4.72 
bees 
7.0) 


1.5 
2),.00 
31.37 


EXPON ( A= 1.75) 


A(t) 


EXPON ( A 


A(t) 


3.51 
4.71 
5.77 
6.7) 


12a0 
1 (eee 
20.51 


(continued ) 


C(t) 


-3h5 
376 
- 105 
139 


2h 
372 


C(t) 


on 
381 
2399 
110 
Lh 
SS10% 
. 230 


C(t) 


-326 
370 
ae lg 
- 386 
337 
21h 
-160 


le? DD 
ELS) 
28.5; 
33.69 
59.8) 
8.09 
95.3) 


21.89 
Bae) 
Wi.1h 
48.59 
82.8) 
97 8h 
SNS 8.8) 





B. DISCUSSION 

The actual probability of crossing appeared to be consistently smaller 
than the theoretical probability under the null hypothesis (A= 1). For 
N of fifty, and over, the agreement appeared to be reasonably good. More- 
over, the probability of rejecting the null hypothesis (power) was quite 
high when A= 1.5 vec was at least twenty-five. This would not be an 
unreasonable sample size for the transistor testing illustration, but 
would undoubtedly by large if the procedure were applied to check entire 
systems. 

Notice that an empirically generated correction could be derived to 
adjust the size of the null hypothesis probability. Consider, for ex- 
ample, the N=h, A= 1.25 element of the nominal five per cent table: 
the crossing per cent (estimated true probability of crossing) is nearly 
five per cent. Under the null hypothesis, the simulation crossing prob- 
ability was computed at 2.3 per cent. Now, if the null hypothesis pre- 
vailed, and if each failure time were multiplied by 1.25)" = 0.8, then 
five percent boundary crossings would be realized. On the other hand, 
if A>1, and each failure time was multiplied by 0.8, but still applied 
against the theoretical derived boundary for A ='1, the power of the 
procedure would be further enhanced. This empirical adjustment has been 
left for more extensive development, e.g. -- for different sample sizes-- 
in the future. It seems clear that an alternate approach involves an 


adjustment of the boundary, itself. 


C. EXPERIMENT 2. 

Using the same procedure as in Experiment 1, we explore the effect 
of «x 4 Bp , but with “PB = constant. We tried the following four com- 
binations: (1) %=1.5B, (2) B=1.5« , (3) X= 2B , (h) B= 2a. 
The results are summarized in the next tables. 
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SABLE iIdal 


5% Probability of a Crossing 
A= 1.9893 B= .99929 


EXPON ( A =.5) 


NC A(t) C(t) f 
B 1.00 W56 Hie, 
1 6.00 1:76 05 
3 (3s 56 Ss 
0 0.00 .000 .00 
0 0.00 .000 00 
0 0.00 .000 00. 
0 0.00 .000 .00 
EXPON ( A= 1) 
NC A(t) c(t) i 
a ). 00 .297 85 
36 5.61 «61 1.80 
bhi 6.95 92 p) OX6) 
ho Be28 50h, 2.45 
hg 17.28 .589 pm li 5 
70 32.91 soo7 3250 
91 7.6 706 55 
EXPON ( A= 1.25) 
NC ~ A(t) C(t) f 
36 4.00 309 1.80 
97 5.59 57 1.85 
118 7.00 93 5390 
17 8.3 2529 7.35 
3h0 18.20 668 16.99 
636 32.99 697 31.79 
871 h6.5h 692 13.54 





= 


NC 


76 
161 
232 
32h, 
857 

1,72 
1799 


NC 


126 
25h 
01 
560 
1381 
1897 
1990 


NC 


169 
02 
59h 
787 
1733 
1987 
2000 


TABLE IIIT (continued) 


EXPON ( A= 1.5) 


A(t) 


4.00 


EXPON ( >= 1.75) 


A(t) 


1.00 
5.55 
6.91 
G27 
to. 6 / 
Doe 
B2a5y 


EXPON ( A 


A(t) 


4.00 
5.50 
6.85 
rar 
15.70 
22.91 
27.50 


C(t) 


316 
28 
495 
ae 1 
.603 
562 
2521 


C(t) 


Bena 
is 


511 
-538 
436 
33h 


C(t) 


mes 
09 


81 
62 
0312 
sae 


al 


3.80 
8.0 
11.60 
16.19 
he 8h 
73-59 
89.9) 


= 


S30 
Jace 
20.0) 
es 
69.0 
9.8) 
99.9 


9 hh 
20.09 
29.70 
39.3h 
86.6), 
993k 


100.00 








TABLE IV 


5% Probability of a Crossing 
A= .99929 P = 1.19893 


= 


EXPON ( A= Sy, 
NC A(t) Cas fy 
10 B.1e ~203 0.50 
hy h.50 Les O20 
10 ane Best 0.50 
3 5.66 ~2)0 0.15 
@) 0.00 .000 0.00 
O 0.00 .000 0.00 
O 0.00 ~000 0.00 
EXPON 

NC A(t) Ce i 

80 3.30 fein ).00 

99 eying 301 he95 

90 Sell 291 1.50 
101 Suey 230); 5.05 

91 12.68 331 1.55 

99 23.06 SS Ton 95 
147 an. 113 ain 5.85 

EXPON ( A = 1.25) 

NC A(t) c(t) ay 
162 3.38 . 263 8.10 
198 Wage 301 9.90 
li 561 323 10.70 
21)2 6.57 Re sell 12510 
392 Shas: Ao ete 19.59 
628 2h. 36 2h13 Bylo, 
oT 35.16 Arey 39.8 
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NC 


21,8 
334 
4,00 
455 


11,06 
Tao 


NC 


346 
477 
588 
(eS 
1368 
1867 
Le ied 


NC 


166 
6h), 
806 
Titel 
171) 
1980 


not run 
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TABLE IV (continued) 


EXPON ( A= 1.5) 
A(t) C(t) 
eee: 268 
heli? 2291 
5.59 sat 
6.58 “358 
13ert 2396 
Oaei)y 386 
31.2h -350 
EXPON ( A = 1.75) 
A(t) C(t) 
3.3 26), 
oly [267 
5.55 Boi 
6260 335 
13.16 361 
20.58 » 303 
Palette 25 
EXPONsiieen = 2 ) 
A(t) C(t) 
3.5 268 
ies 293 
Sos 307 
6.53 323 
12eeT 316 
17.39 216 


12.39 
16.70 
Se 
22.7h 


70.29 
85.99 


Vises 
23.8), 
29.39 
36.15 
68.39 
93 3h 
98.95 


23-29 
32.19 
0.29 
18.8), 
85.70 
98.99 








TABLE V 


5% Probability of a Crossing 
X= 1.73082 B= .865))1 


EXPON ( A= .5) 
NC A(t) C(t) fi 
2 4.00 .156 6x10 
1 6.00 0:76 0.05 
1 8.00 6119 0.05 
0 0.00 .000 0.00 
O 0.00 000 0.00 
0 0.00 .000 0.00 
0 0.00 .000 0.00 
EXPON ( A= 
= A(t) CG) t 
h 00 .169 0.20 
17 5.88 hhs 0.85 
28 7.50 56h) 1.0 
3h sie ~507 1G 
39 118,92 702 1.95 
59 36.86 858 2.95 
85 52.07 8S y.25 
EXPON ( A= 1.25) 
NC A(t) C(t) fy 
11 1.00 -207 0.55 
18 5.87 1460 2.h0 
70 Peo 538 3.50 
103 8.8) 568 5.15 
308 19.8 ~805 15.h0 
612 35.67 .808 30.59 
851 50.50 803 2 .5h 
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TABLE V (continued) 


EXPON ( )\= 1.5) 


{=a 


NC A(t) CG) i 
any 1.00 .1958 0.85 
9h 5.86 ohh its 168; 
161 Ae, 53h 8.0; 
253 8.89 pees 12.6); 
80), 18.82 691 10.19 
1163 327, );2 651 731k 
1802 3.71 .589 90.09 

EXPON ( A = 1.75) 

NC A(t) OC) fy 
28 h.00 20) 1.0 
169 5.81 sey 8.5 
286 7B 8 fay 1.29 
W)i7 8.83 571 23h 
1327 18.03 .620 66.3), 
1895 28.73 192 9. 7h 
1987 Bone 381 99.3) 

EXPON ( d= 2) 

Nc A(t) CG) fy 
Sh 00 MPG 2 
239 5.81 18 11.95 
M7 fieele «97 22.3h 
666 Gaf0 538 33.29 
en A 2 2535 86.20 
1985 25.28 - 358 99.2h 
2000 30.7h .265 100.00 
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TABLE VI 


5% Probability of a Crossing 
a= .865)1 B= 1.73082 


NC A(t) C(t) fy 
1h 2.78 Gil D6 
6 4.16 258 0.30 
i) 5.5 31h 0.55 
5 80 .176 0.25 
O 0.00 .000 0.00 
1 3.00 .130 0.05 
1 8.00 .130 0.05 

EXPON ( A= 1) 

Nc A(t) CCE) -, 

Hk 3.02 5208 1.85 
Tae? 3.95 e230 5.85 
106 5.00 225) 5.30 
late 5.60 sa Bohs 
105 11.48 261 5.25 
119 20.31 .30) 5.95 
122 30.05 336 6.10 

EXPON ( A= 1.25) 

NC A(t) C(t) if 
192 3.10 (2ec 9.60 
228 Inaile 2253 11.0 
25 5.10 27h eel 
265 5.98 »279 13205 
1,00 12312 328 19.99 
600 Zine © 337 la ENS) 
(77 So ale ans Wh. 38.8) 


2h 








NC 


299 
361 
22 
81 
859 
1363 
1679 


NC 


526 
682 
837 
991 
1703 
1976 


not run 


TABLE VI (continued) 


EXPON ( »= 1.5) 


BXPON 


tr 


s— 
ig 


ati 
=) 
me 


25 


(We ASS) 


= PP) 


C(t) 


232 
+253 
269 
287 
323 
3) 
312 


Gey) 


226 
25h 
-263 
281 
. 309 
271 
212 


c(t) 


0229 
250 
262 
280 
2278 
19) 


1h.95 
18.0) 
21.09 
2),.0h 
42.95 
68.1); 
83.94 


19.9) 
25.20 
30.6) 
36.39 
65.99 
92.1) 
98.59 


26.29 
3h.09 
heel 
ho .5h 
85.1) 
98.79 








D. DISCUSSION 

The tendency of the theory to understate the crossing probability, 
noted “= $ , is even more pronounced when 4= 1.58. However, when 

B=1.5% (TABLE IV) the understatement error nearly vanishes, and the 

test is very close to being exact (i.e.- crossings occur with probability 
about five percent) oe all values of N checked. According to TABLE VI, 
the probability of rejection is slightly overstated when B= 2%, for 
all but N = ). To within limits of sampling accuracy, it appears that 
PB equalling 1.65 will give a very nearly correct answer at 4 = 0.05. 
Further sampling will pin this down more precisely. A theoretical ex- 
planation would be most interesting, for the asymptotic (large N) mathe- 
matics detects no difference in situations for which A f= constant. 

The next set of experiments was conducted in order to investigate 
our procedure's ability to detect departures from the null (exponential, 


A= 1) hypothesis. 


BK. EXPERIMENT 3. 

Exvloration of the crossing probability when the boundary reflected 
the null hypothesis (exponential with A = en cee “<=B8 , nominal 
crossing probability five percent), but samples were drawn from a special 


long-tailed distribution with a Pareto-like appearance: 


ax 





B(x) = ee 0. (4.1) 


|] +ax 


In order to give a fair comparison, the parameter a was chosen so as to 
match the medians of i and the null distribution; thus a = Ge)e 
A rationale for the choice of (h.1) is that such a distribution 
appears very Similar to the exponential, but has a decreasing failure 

rate and a bias toward early failures. The results are recorded in 


TABLE VII. 
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TABLE VIL 


5% Probability of a Crossing 
X= B = 1.22388 


. PARETO 

N NC A(t) C(t) ee 
h 19 3mibO 1253 3.95 
6 97 ).82 2319 h.85 
8 91 5.61 Se h.55 
10 10h) 6.70 ~ 309 5.20 
25 133 12.30 eer 6.6), 
50 209 20.00 257 10.5 
75 218 26.65 aie 10.90 


F. DISCUSSION. 

Our sampling experiment indicates that the sequential procedure under 
study has very little power against the "Pareto tail" alternative ().1) 
under the particular choice of parameters, at least for small N. Actually, 
the probability of rejection is somewhat larger than for the correspond- 
ing exponential. It seems likely that the relative prevalence of small 
observations accounts for the increased crossing probability, and also 
for the relatively short expected time to cross, given that a crossing 


occurred. 


G. EXPERIMENT l,. 
The procedure used in Experiment 3 was repeated, except that samples 


were drawn from a shorter-tailed (Weibull, with increasing hazard rate) 


distribution: 
ee 


F(x) = 1-e ; (.2) 


Again, the parameter a was fitted by matching medians of the exponential 


val 





( >° ame => 


—— « 





~1 
to (4.2); by coincidence, a = (log,2) once again-- an accident of 


parameterization. The results are recorded in TABLE VIII. 


TABLE VIILI 
5% Probability of a Crossing 
a= B= 1.22388 
N NC A(t) C(t) i 
h 13 leee 25 0.65 
6 38 5.89 .610 ace 
8 58 Vane - 702 2.90 
10 91 9.6 ate 4.55 
25 805 2270 1.078 b@nen 
50 189), hi .88 1.091 9.69 
15 1996 58.30 1.009 99.79 


H. DISCUSSION. 

The sampling experiment suggests that the power of our sequential 
procedure is quite unsatisfactory against the particular Weibull alterna- 
tive examined for small N, but the power rapidly increases with N. Scru- 
fmaizine of the density of ()).2) indicates that there are relatively few 
"early failures", while a concentration of "wearouts" probably accounts 
for the rather long conditional expectation of the time to hit the 


boundary (reject the null hypothesis). 
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